Novel and simple description for smooth transition from a cluster to jj-coupling shell 

model wave function 



o 

(N 
X> 

IX, 
m 

(N 



Tadahiro Suhara 1 , Naoyuki Itagaki 2 , Jozsef Cseh 3 , and Marek Ploszajczak 4 
1 Institute of Physics, University of Tsukuba, Tsukuba 305-8571, Japan 
2 Yukawa Institute for Theoretical Physics, Kyoto University, 
Kitashirakawa Oiwake-Cho, Kyoto 606-8502, Japan 
3 Institute of Nuclear Research of the Hungarian Academy of Sciences, Debrecen, Pf. 51, Hungary-4001 
4 Grand Accelerateur National d'lons Lourds (GANIL), 
CEA/DSM - CNRS/IN2P3, BP 5027, F-14076 Caen Cedex 05, France 
(Dated: February 26, 2013) 

We propose an improved version of Antisymmetrized Quasi-Cluster Model (AQCM) to describe a 
smooth transition from the a-cluster wave function to the jj'-coupling shell model wave function and 
apply it to the ground state of 12 C. The cluster-shell transition in 12 C is characterized in AQCM by 
only two parameters: R representing the distance between a clusters and the center of mass, and A 
describing the break of a clusters. The optimal AQCM wave function for the ground state of 12 C 
is an intermediate state between the three-a cluster state and the shell model state with the P3/2 
subshell closure configuration. The result is consistent with that of the Antisymmetrized Molecular 
Dynamics (AMD), and the optimal AQCM wave function quantitatively agrees with the AMD one, 
although the number of degrees of freedom in AQCM is significantly fewer. 

PACS numbers: 21.30.Fe, 21.60.Cs, 27.20.+n, 21.60Gx 
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I. INTRODUCTION 

Shell structure is a fundamental property of atomic nu- 
clei. The stability of nuclei and presence of magic gaps is 
explained by non-uniformities of the single-particle level 
distribution The strong spin-orbit interaction is es- 
sential to explain the observed magic numbers 0, Q ■ The 
nuclear shell model in which the one-body behavior is 
supplemented by configuration-mixing effects of residual 
two-body interaction, describes nucleus as a closed quan- 
tum system where nucleons occupying bound orbits are 
isolated from the environment of scattering states and 
decay channels. In its modern version, shell model calcu- 
lates nuclear properties in ab initio manner Q , using re- 
alistic interactions which reproduce the nucleon-nucleon 
scattering data [H, Q . 

The validity of such a closed quantum system 
framework depends on the dissociation energy. The 
configuration-mixing effects in weakly bound or unbound 
nuclear states cannot be treated as a small perturba- 
tion atop the mean field, and involve effects of the cou- 
pling to decay channels 0j- In particular, the appear- 
ance of cluster states in the vicinity of their respective 
cluster decay thresholds is a consequence of an openness 
of the nuclear many-body system [H, . Consequently, 
the standard shell model approach simply cannot be suc- 
cessful in predicting cluster states at low excitation en- 
ergies around the cluster-decay thresholds. In the most 
advanced closed quantum system approaches to cluster 
decay, the shell model wave functions must be supple- 
mented with a cluster component to achieve a quantita- 
tive agreement with the data fioj ). 

The failure of the closed quantum system approaches 
to describe cluster states is the central problem in nuclear 
theory. The energetic order of particle emission thresh- 



olds, and their nature, depends on precise properties of 
the nuclear Hamiltonian. On the other hand, the phe- 
nomenological rule that cluster correlations are seen only 
in the vicinity of the respective cluster emission threshold 
cannot be a consequence of specific properties of nuclear 
forces. A generic explanation of this rule in terms of the 
collective external coupling of shell model states via the 
decay channel(s) has been put forward in Refs. 0, 
In this context, a-cluster states are of particular interest 
because of strong binding of an a particle and a weak a-a 
interaction which does not allow to bind an a-a system. 
The systematics of a-cluster states in light nuclei was a 
basis for Ikeda conjecture that a-cluster states appear 
close to their cluster decay thresholds . 

In spite of recent remarkable advances in the open 
quantum system formulation of the shell model [T2I [l3j j , 
a coherent picture of shell structure and clustering in the 
framework is not yet within our grasp. The closed quan- 
tum system description of low-energy excitations in light 
nuclei, where cluster- and shell-model-structures are in- 
tertwined, requires a generalization of the shell model 
wave function. 

On the other hand, cluster models have been devel- 
oped and successfully applied to describe cluster states. 
In traditional cluster models, each a cluster is expressed 
as (si/2) 4 shell model configuration. In this case, the 
spin-orbit interaction, which plays an essential role in 
the nuclear systems, cannot be taken into account. This 
failure is coming from the special symmetry of the wave 
function. Since four nucleons shares the same spatial 
wave function, each a cluster is a spin-zero system, and 
the spin-orbit interaction does not work. If we take 
the limit of zero distance among the clusters, the clus- 
ter model wave function agrees with the shell model one 
in the case of closed shell nuclei such as 16 O and 40 Ca. 
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However, it cannot describe a subshell closure configura- 
tion, where the contribution of the spin-obit interaction 
becomes maximum. Therefore, extension of the clus- 
ter model space, especially for the spin configurations, 
is needed for the general description of the nuclear struc- 
ture. In recent years, there have been many microscopic 
studies along this line, and competition of cluster and 
shell components in the ground state of light nuclei have 
investigated [2 EH- For example, the ground state of 
12 C is an intermediate state between the three-a cluster 
state and the shell model state with the p 3 / 2 subshell 
closed configuration However it is still dif- 

ficult to evaluate quantitatively to what extent cluster 
structures develop or shell model structures admix. 

A step forward in this direction is proposed in the 
Antisymmctrizcd Quasi-Clustcr Model (AQCM) (lii - E^ 
which attempts to include these two distinct structures of 
different physical origins in a single many-body approach. 
In the AQCM, the transition from the cluster- to shell- 
model-structure can be described by only two parame- 
ters: R representing the distance between a clusters and 
the center of mass, and A which characterizes the quasi- 
cluster(s) and quantifies the role of the spin-orbit interac- 
tion in breaking the a cluster (s). This is very transparent 
and simple approach to quantitatively discuss the mixing 
of cluster and shell components. However, the previous 
AQCM, which were applied to the cluster-shell compe- 
tition in C, Ne and Mg isotopes [l9l - [23j . has a problem; 
description for the subshell closure configurations was not 
exact. For example, in the studies of C isotopes (l9l. [20|. 
only one of the three a clusters was changed into a quasi- 
a cluster which corresponds to the P3/2,±3/2 shell model 
orbits. The P3/2,±i/2 shell model orbits were not included 
in the model space because the remaining a clusters were 
unchanged. The purpose of the present work is to im- 
prove the AQCM description by removing the restriction 
for the spin orientations of individual nucleons. In the 
new formulation, all a clusters can be changed into quasi- 
a clusters, and jj-coupling shell model states including 
the subshell closure configuration can be described. 

The paper is organized as follows. The formulation 
and details of the improved AQCM parameterization are 
given in Sect. |H] In Appendix [XJ we prove that this 
AQCM can describe the p$/2 subshell closure configura- 
tion. In Sect. IIIH the AQCM results for 12 C arc discussed 
and compared with the results of the Antisymmctrizcd 
Molecular Dynamics (AMD). We also discuss the rela- 
tion of the AQCM with the SU(3) shell model. Finally, 
in Sect. IIVI we summarize the results and give the main 
conclusions. 



II. THE MODEL 

In this section, wc discuss the improved parameteri- 
zation of the AQCM wave function and the many-body 
Hamiltonian used in this work. 



A. Single particle wave function (Brink model) 

In conventional a cluster models, the single particle 
wave function is described as a Gaussian packet [24| : 
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where rji represents the spin-isospin part of the wave func- 
tion, and _R 7 is a real parameter representing the center 
of a Gaussian for the 7-th a cluster. In this Brink-Bloch 
wave function, four nucleons in the 7-th a cluster share 
the common R 1 value. Hence, the contribution of the 
spin-orbit interaction vanishes. 



B. Single particle wave function in the AQCM 

In the AQCM, a nucleus consists of quasi-a clusters. 
For nucleons in the quasi-a cluster, the single particle 
wave function is described by a Gaussian packet, in the 
same way as in the Brink-Bloch wave function. However, 
the center of this packet is a complex parameter 
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Xi and Ti in Eq. @ represent the spin and isospin part of 
the i-th single particle wave function, respectively. For 
the width parameter, we use the value v = 0.235 fm -2 . 
The spin orientation is given by the parameter while 
the isospin part is fixed to be 'up' (proton) or 'down' 
(neutron), 



* = &tlt)+£uU>. 

Ti = \p) or 1 77.). 



(4) 
(5) 



In Eq. ([3]), e^ pin is a unit vector for the intrinsic-spin 
orientation, and A is a real control parameter describing 
the dissolution of the (quasi)-a cluster. As one can see 
immediately, the A = AQCM wave function, which has 
no imaginary part, is the same as the conventional Brink- 
Bloch wave function. We explain later that the AQCM 
wave function corresponds to the j j-coupling shell model 
wave function when A = I and i2 7 — ¥ 0, and the im- 
proved AQCM approach can describe the subshell closure 
configuration. 

The spin-orbit interaction is intuitively interpreted as 



I ■ s = (r x p) ■ s = (s x r) • p 



(0) 



where r, p, and s are the position, the momentum, and 
the spin of the nuclcon, respectively. If nucleons have the 
momentum components parallel (anti-parallel) to s x r, 
the spin-orbit interaction acts attractively (repulsively). 
An imaginary part of a Gaussian wave packet means the 
momentum of nuclcon: 



{A\Pi\i>i) = 2^v^Im[0 



(7) 
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Therefore, the AQCM wave function has the momentum 
parallel or anti-parallel to s x r by introducing this partic- 
ular form (Ae^ pm x _R 7 ) of the imaginary part in Eq. (|3j). 
A > and A < correspond to parallel and anti-parallel 
momenta, respectively. In actual calculation, of course, 
the spin-orbit interaction is a two-body force, which is a 
function of — rj and not of r. 

In this paper we focus on 12 C, which consists of three 
a clusters, and we assume that they are placed with an 
equilateral triangular configuration. For the first quasi- 
a-cluster, the spin direction is defined along the z-axis 
and the center of mass is set in the x direction as Re x . 
Then we introduce imaginary parts of centers of Gaus- 
sians related to the momenta of nucleons as in Eq. ([3]), 
which means in the y direction. Hence, the centers of 
Gaussian wave packets are 

C = R{e x + ihey) (8) 

and 

d = R(e x - iAe y ) (9) 

for the spin-up and spin-down nucleons, respectively. 
Here, e x and e y are unit vectors on x and y axes, re- 
spectively. 

For the second and third quasi-a-clusters, we construct 
their spin directions and the centers of Gaussian wave 
packets by rotating the first quasi-a cluster. Since both 
the spatial and spin parts of the wave function are rotated 
simultaneously, the relative angles among r, p, and s in 
the first quasi-a-cluster are kept in the second and third 
quasi-a clusters. Thus, the spin-orbit interaction also 
acts in these quasi-a clusters as well as in the first one. 
We rotate both the spatial and spin parts of the first 
quasi-a-cluster about the y-axis as 

V>i+4 = £(<* = 0,0 = 01,7 = 0)^, (10) 

^4+8=^(0 = 0,0 = ^2,7 = 0)^, (11) 

where i — 1 — 4, (a, 0,7) = f2 are the Euler angles, 
R(fl) = e~ laJ *e~ l P Jy e~ l ' yJ z is the rotation operator, and 
6\ and 6*2 are rotational angles, 6\ = 2tt/3 and 6* 2 = 
47r/3. From the discussion in the next subsection, it is 
understood that the AQCM wave function coincides with 
the (si/2) (P3/2) 8 configuration under the proper 6±, 62 
values and the conditions of A = 1 and R — > 0. 

The essential difference between the previous and 
present AQCM comes from the treatment of the spin ori- 
entation. In the early version of the AQCM [19l - l21j |. the 
intrinsic spin of each nucleon was quantized with respect 
to the z-axis and the spin direction was restricted to z or 
— z direction for all nucleons 

i= (1,0) or (0,1) (12) 

Therefore, it was impossible to change all a clusters to 
quasi-a clusters and describe the subshell closure config- 
uration. 



In the improved AQCM, we expand the model space 
and spin directions can be oriented in any direction. This 
means £^ and £4 are complex parameters depending on 
the rotational angles Q\ and 62- As a result, we can 
change all a clusters to quasi-a clusters by rotating both 
the spatial and spin parts of the wave function of the 
first quasi-a-cluster, and describe the P3/2 subshell clo- 
sure configuration. 

1- P3/2 subshell closure configuration 

Let us now discuss the AQCM single particle wave 
function analytically to prove that it describes the P3/2 
subshell closure configuration at R — > and A = 1 . 

First, we discuss the spin- up nucleons in the first quasi- 
a cluster, whose centers of Gaussian wave packets are 
given in ([5]). The single particle wave function is: 
3 

^ ^ ( 4 cx p {- vr l - < 2 + 2vr * ■ Ci] 1 1 )n. (13) 

Using Eq. ©, the last factor can be expanded as 
cxp [2vn ■ C] = fcj (2^(*i + «AyO) fc , 

k=0 

For A = 1, using the relations 

/^+%\ fe = L Ykk (n i ), (is) 

\ n J s k 

and 

rf cxpf-^r 2 ] = j-Rokin), (16) 
the single particle wave function becomes 

^ = (—) 4 E §^-Rofc(r<)n*(fii)| t )n. (17) 

Here, s k and t k are the normalization factors of spherical 
harmonics Ykk(^i) and the radial wave function Rok(fi), 
respectively. Since the spin is up, this wave function 
is described as a linear combination of \j,j z = j) = 
Rok{r l )Y kk (n i )\ | ) states where j = k + 1/2: 

OO 

&= a i Rj ~ h \iii)Ti- (is) 

J = l/2 

a,j in (I18p is the coefficient for the \j,j) state with the 
separated factor of R^^ . Analogously, the single parti- 
cle wave function for the spin-down nucleons in the first 
quasi-a cluster can be written as 

00 

V^= a-jR'-^-iW (19) 

J=l/2 
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This expression can be obtained by a time-reversal trans- 
formation of the spin-up wave function. 

Next, we discuss nucleons in the second and third 
quasi-a-clusters which are generated by multiplying the 
rotational operator for the first quasi-a cluster. Accord- 
ing to the definition of Wigner D function, 

(j,k\R(a,P,i)\j,m) = exp[-ika}d 3 km (f3) expl-imj] 

(20) 

where d 3 km {(3) is a Wigner small d function, the rotated 
state can be expressed as 



£(a = 0,0 = 0,7 = 0)|.7\ j) 



j 

E 



d^(0)b\m)- (21) 



Thus, the single particle wave functions in the second 
and third quasi clusters, which are generated by rotating 
\j,j) in the first quasi-cluster about the y-axis, can be 
written as 

£(a = O,/3 = 0,7 = O)Vi 

oo j 

= E E ajR'-l^OMmfr. (22) 

j=l/2 m=-J 

We perform this transformation also for the spin-down 
nucleons in the first quasi-a cluster, which are expressed 
as linear combinations of {\j,—j)}- The rotation mixes 
other \ j, m) states than the \j, ±j) states, when satisfies 

Rotating all single particle wave functions following 
Eq. (f2"2"|) and substituting them into Eq. (|23p. we can 
discuss the nature of the Slater determinant consisting 
of quasi-clusters. Owing to the anti-symmetrization, the 
lowest order of the wave function in R is R 8 (R 4 for 
both proton and neutron parts), and the coefficient for 
R 8 exactly corresponds to the (si/2) 4 (.P3/2) 8 configura- 
tion. Therefore, for A = 1, R — > 0, and proper 9\ and 02 
values, the total wave function coincides with the Ps/2 
subshell closure configuration (see Appendix IS]) . This 
limit is different from the R — > limit of the three-a 
cluster model, which is known as the Elliott SU(3) limit, 

( S ) 4 ( Px ) 4 (Py) 4 . 



C. AQCM wave function of the total system 

The AQCM wave function of the total system consist- 
ing of three quasi-a-clusters is projected onto parity and 
angular momentum eigenstates using the parity projec- 
tion operator and the angular momentum projection 
operator P( [K as 

® = Y, CkP mk P± * / ^^--- >^ia] (23) 

K 

The parity projection operator P ± is defined as 



1±P 



where P is the parity operator. The angular-momentum 
projection operator P^k ^ s defined as 



pJ 



2J + 1 
8tt 2 



dnD J A i K {n)R(n), 



(25) 



where Df IK (fl) is the Wigner D function. In the total 
wave function 



Eq. (|23|) . particles whose numbers are 
from AN — 3 to AN form the iV-th quasi-a cluster. 



D. Hamiltonian 

The Hamiltonian operator H has the following form: 

A A 



H = y~ y Jj- T c . m . + ^ Vij > 



(26) 



i>J 



where a two-body interaction &j includes the central 
part, the spin-orbit part and the Coulomb part. For the 
central part, we use the Volkov No. 2 effective N — N 
potential |25| : 



^ c (^) = E^ ex p(- 



rt 3 /ci)(W 



MP a P T ), (27) 



fe=i 



where M = 0.60, W = 1 — M = 0.40, = -60.65 MeV, 



Vi 



61.14 MeV, a = 1.8 fm, and c 2 = 1.01 fm. For 



the spin-orbit potential, we adopted the spin-orbit part 
of the G3RS potential (26|: 



V 1 



E^ S 
fc=i 



cxp(-f?./d 2 fe ) J P( 3 0) J L • S, (28) 



(24) 



where V( s = -1600 MeV, V\ s = 1600 MeV, d x = 0.600 
fm, c?2 = 0.477 fm, and P( 3 0) is a projection operator 
onto a triplet-odd state. The operator L stands for the 
relative angular momentum and S is the spin operator. 

The parameter set: M = 0.60, V-} s = -2000 MeV, and 
V\ s = 2000 MeV is known to give a reasonable descrip- 
tion of a + 7i and a + a scattering phase shifts ]27| . In 
the present study, we adopt, however, the slightly weaker 
strengths of the spin-orbit force to fit the 0^~ energy of 

l2 c m. 



III. RESULTS AND DISCUSSION 

In this section, we show results of the AQCM for 
12 C and compare them with the AMD results. In 
particular, we extract the single particle orbits with 
the Antysymmetrized Quasi-Cluster + Hartree-Fock 
(AQC+HF) method to show that the improved AQCM 
can describe the P3/2 subshell closure configuration. We 
also discuss the relation of the AQCM with the SU(3) 
shell model. 
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FIG. 1: (Color online) The energy surfaces for the + state 
of 12 C, the total energy (top) and the spin-orbit energy (bot- 
tom), respectively. 

A. Energy surfaces 

The AQCM (three quasi clusters) energy surfaces of 
12 C as functions of R and A are shown in Fig. [TJ The 
top and bottom panels show the energy and spin-orbit 
energy surfaces for the + state of 12 C, respectively. 
The minimum of the energy surface appears at around 
(R, A) = (0.9, 0.2) and the corresponding energy is —89.6 
MeV. This small distance between clusters, R = 0.9 fm, 
and a small but finite imaginary part of the wave func- 
tion, A = 0.2, indicate that the ground state of the 12 C is 
an intermediate state between the three-a state (A = 0) 
and the p 3 / 2 subshell closure state (R — s- and A = 1). 
This hybrid character of the ground state wave function 
can be confirmed by calculating the squared overlap of 
the wave functions: 

A) = (0.9,0.0))|$ min )| 2 = 63.6%, (29) 
|($((i?,A) = (0.01, 1.0))|$ min )| 2 = 47.3%, (30) 

where $ m i n is the wave function at the minimum of the 
energy surface. One can see that the squared overlaps 
between $ m ; n and the three-a state ((R,A) = (0.9,0.0)), 
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FIG. 2: The AQC+HF effective single particle energies of the 
improved AQCM are plotted as a function of the parameter A 
(R = 0.01 [fm]) in 12 C. The Coulomb interaction is neglected 
in this calculation. 



and between $ m i n and the P3/2 subshell closure state 
((R, A) = (0.01, 1.0)) have both significant values around 
0.5. 

From the spin-orbit energy surface, one finds that the 
spin-orbit interaction acts attractively and repulsively for 
A > and A < 0, respectively (for A = 0, the spin-orbit 
interaction does not act). Around the minimum of the 
total energy, the spin-orbit interaction contributes ^ — 16 
MeV to the total energy. In the next subsection, we 
will show that the improved AQCM takes into account 
effects of the spin-orbit interaction very well, and the 
AQCM wave function is almost the same as the AMD 
wave function. 



B. Single particle orbits of the AQC+HF method 

The AQCM single particle wave functions are given 
by Gaussian packets with complex parameters, and dif- 
ferent single particle wave functions arc mutually non- 
orthogonal. To extract orthogonal single particle or- 
bits and corresponding effective single particle energies 
from AQCM, we apply the method imitating Hartree- 
Fock (HF) approach which was applied first in the con- 
text of the AMD 0. We call it the AQC+HF method. 
In general, the corresponding AQC+HF single particle 
orbits are given by linear combinations of Gaussian sin- 
gle particle wave functions. One should stress however 
that the HF self-consistency in AQC+HF approach is 
satisfied only within the restricted functional space of 
AQCM single particle wave functions. Moreover, except 
for special cases, resulting single particle orbits are not 
eigenstates of angular momentum and parity operator. In 
spite of these restrictions, the AQC+HF approach pro- 
vides useful information about the dependence of mutu- 
ally orthogonal single-particle orbits and effective single 
particle energies on parameters of the AQCM manifold. 
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FIG. 3: The same as in Fig. [2] but for the old version of the 
AQCM [HQ. 



In the following, we analyze these effective single particle 
energies to show that the improved AQCM describes the 
P3 li subshell closure with the set of the parameter values 
R^0 and A = 1. 

Fig. [5] shows AQC+HF effective single particle ener- 
gies of 12 C for the improved AQCM as a function of A. 
The Coulomb interaction is neglected and the limit of 
R — > is expressed by R = 0.01 fm. There are in this 
case only three AQC+HF single particle orbits expressed 
by solid, dashed, and dotted lines. Each of these orbits 
is occupied by four nucleons because of the spin-isospin 
degeneracy. The lowest orbit (the solid line) corresponds 
to S1/2 shell model wave function. The higher two or- 
bits (dashed and dotted lines) are linear combinations 
of p 3 /2 and P1/2 components. At A = 1, these two sin- 
gle particle orbits become degenerate. In the vicinity of 
A = 1, the dominant component in higher orbits is Ps/2- 
The P1/2 component vanishes at A = 1. At this point, 
the higher single particle orbits are pure p^/2 shell model 
wave functions. If we take into account the spin-isospin 
degeneracy of four, eight nucleons occupy these degen- 
erate ps/2 orbits. Hence, for R — > and A = 1, the 
AQCM wave function describes the ground state con- 
figuration (si/2) 4 (P3/2) 8 of 12 C in the jj-coupling shell 
model. For other A values, there are three independent 
AQC+HF single particle orbits. In this case, the AQCM 
wave function does not describe the p 3 / 2 subshell closure 
configuration. 

Fig. [3] shows the AQC+HF effective sin gle particle 
of 12 C for the old version of the AQCM where 
only one a cluster is changed into a quasi-a-cluster. In 
this case, there exist always three independent AQC+HF 
single particle orbits for any value of A, i.e. the old 
version of the AQCM cannot describe the P3/2 subshell 
closure configuration. 

In Fig. 21 we show the expectation value of the one- 
body spin-orbit operator, Li ■ Si, as a function of R 
and A. Here the sum over i is for all nucleons. The cal- 
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FIG. 4: (Color online) The expectation value of one-body 
spin-orbit operator, J2i Li ■ Si, in the improved AQCM. 

culations are performed using the improved AQCM. The 
expectation value becomes zero for the three-a model 
(A = 0), and its sign changes when crossing the A = 
line. For the P3/2 subshell closure configuration, the 
expectation value of the one-body spin-orbit operator 
should be 4, i.e. 8 (number of nucleons in the P3/2 or- 
bit) x 1/2 (the expectation value of Li ■ Si for the p 3 / 2 
orbit). Indeed, we can see that this value is reached ex- 
actly in the limit of A = 1 and R — >• 0. This proves again 
that the AQCM can describe the P3/2 subshell closure 
configuration of the jj-coupling shell model. 



C. Comparison of AQCM and AMD results 

In this subsection, we compare our results of the 
AQCM with AMD. In Rcf. [H, we have calculated 12 C 
using the f3-j constrained AMD with the same interac- 
tion as used in the present work. In AMD, all nucleons 
are described as independent Gaussian wave packets, and 
the Gaussian center parameters are complex variational 
parameters. This variation called 'cooling process' is of- 
ten performed before the angular momentum projection. 
Here, we have introduced as constraints the quadrupolc 
deformation parameters, (3 and 7, and prepared many dif- 
ferent states to solve the cooling equation. By projecting 
all these different configurations onto the eigenstates of 
angular momentum and parity, we can obtain the lowest 
energy configuration after the projection. Thus the /3-j 
constrained AMD somehow overcomes the approxima- 
tion of projection after variation and can be considered 
as an improved version of AMD. 

We compare the optimal solution of our AQCM wave 
function (R, A) = (0.9,0.2) with this constrained 
AMD at the minimum point of the + energy surface on 
the (9-7 plane, ((3 cos 7, (3 sin 7) = (0.35,0.17), which is 
the dominant component of the Of state of 12 C. The en- 
ergy associated with these AQCM and AMD wave func- 
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tions arc —89.6 McV and —90.1 McV, respectively, i.e. 
their difference is rather small, 0.5 MeV. In addition, 
the squared overlap between the AQCM and AMD wave 
functions is very large, 



|($min|$A 



MD.min/ 



= 98.6%, 



(31) 



where $AMD,min is the + projected AMD wave func- 
tion at (j3 cos 7, (3 sin 7) = (0.35,0.17) in Ref. [H|. This 
very large overlap indicates that the AQCM wave func- 
tion gives an almost identical result of the AMD wave 
function for the ground state of 12 C. Remaining small 
differences between AMD and AQCM wave functions are 
due to different symmetries in both models as discussed 
later. The AQCM wave function has a threefold rota- 
tional symmetry about the j/-axis, while the AMD wave 
function has no symmetry 

The number of degrees of freedom in AQCM and AMD 
wave functions is very different. The AQCM wave func- 
tion has only two degrees of freedom (R and A) , whereas 
the AMD wave function in the studied case has ~ 100 
degrees of freedom, such as real and imaginary parts of 
the center of Gaussian wave packet and the direction of 
the spin for each single particle wave function. This in- 
dicates that the AQCM describes effect of the spin-orbit 
interaction and the breaking of a clusters very efficiently 

In the framework of constrained AMD, wave func- 
tions with the constraint of 7 = 7r/3 have approximately 
threefold rotational symmetry. Therefore, one can com- 
pare AQCM and AMD wave functions having the same 
threefold symmetry by taking the minimum AMD wave 
function on the 7 = 7r/3 line. This minimum point 
is 7) = (0.29, 7r/3) and the corresponding energy is 
—89.6 MeV, which is the same value of the optimal so- 
lution of the AQCM wave function. In addition, the 
squared overlap is very close to unity, 



$AMD )| 2 = 99.9%, 



(32) 



where $AMD, sy m is the + projected AMD wave function 
at (/3, 7 ) = (0.29,7r/3). This almost 100 % overlap of 
these two wave functions means that the AQCM result 
is consistent with the AMD result for the ground state 
of 12 C except for the effect of symmetry breaking. 

Finally, we compare the ,8-7 constrained AMD with 
the previous AQCM wave function [l9l - [2l| . where only 
one of the a clusters is changed into the quasi-cluster. 
In this version of the AQCM, the minimum point of the 
energy surface is at (R, A) = (1.2,0.4) and its energy is 
—89.2 MeV. Then, energy difference between the previ- 
ous AQCM and AMD is 0.9 MeV, which is not a serious 
problem compared with the improved AQCM, 0.5 MeV. 
However, the squared overlap of these two wave functions 
deviates very much from unity 



pre, mm 



*AMD,min)| 2 =90.8%, 



(33) 



where & r 



is the previous AQCM wave function at 



the minimum of the total energy. 



D. Relation with SU(3) shell model 

The relation between the shell and cluster models was 
established historically in the harmonic oscillator limit 
[30| . i.e. via the SU(3) symmetry [SJ [Hj], and much of 
the cluster-shell coexistence or cluster-shell competition 
has been investigated along this line. Our present study 
in the AQCM framework represents a different approach, 
and therefore, it is an interesting question, how the two 
methods and their results compare to each other. 

The harmonic oscillator interaction, and the SU(3) ba- 
sis states represent the intersection of the shell and clus- 
ter models [30l l3l| . Therefore, in the symmetry-adopted 
models, like SU(3) shell model [32|, or its multi-shell ex- 
tension, the symplectic shell model (33| . on one side, or 
the semimicroscopic algebraic cluster model [34[ on the 
other side, it is easy to characterize the relative impor- 
tance of the shell and cluster components. Actually, there 
are two straightforward ways to do so. One of them is to 
determine the SU(3) content of the realistic wave func- 
tion, i.e. to expand it in terms of the SU(3) basis. (The 
shell [H, H|| and cluster [34| models mentioned here ap- 
ply SU(3) basis, thus this is an easy task.) In the extreme 
case of a single SU(3) basis state, the overlap is 100 per- 
cent, and in general the leading representation shows the 
overlap, i.e. how much the state is a shell-like cluster 
state. The other way is to split the interaction into an 
SU(3)-preserving and a symmetry-breaking part. Their 
relative weights can be used as a control parameter to 
indicate how close the real situation is to the intersec- 
tion point of the two models [35| . This latter method is 
applied in the studies of the shape-phase-transitions of 
nuclei (35l. l36j. 

The Hamiltonian can be separated easily for some 
model interactions. In the shell model, typical 
symmetry-preserving parts are the harmonic oscilla- 
tor and quadrupole forces, while important symmetry- 
breaking parts are the spin-orbit and pairing interac- 
tions. In the cluster model, the symmetric parts arc again 
the harmonic oscillator and quadrupole forces, and the 
symmetry-breaking interaction can be e.g. the dipole in- 
teraction. 

In the Op shell, the SUY3) model reduces to Wigner's 
supermultiplct theory 37J, and therefore, it is a very- 
well studied region. The joint conclusion of many works 
is that for the ground state of the 12 C nucleus the (0,4) 
SU(3) symmetry is a good approximation. The recent 
developments using the no-core shell model HH confirms 
the dominance of the (0, 4) component in the wave func- 
tion m. 

The AQCM proposed in this work, provides a pos- 
sibility for going from the molecule-like clusterization 
to the jj coupled shell model limit directly, in a sin- 
gle step, contrary to the two-step procedure of the shell- 
and cluster-model calculations with different symmetry- 
breaking terms, as we discussed in Ref. (22[. This model 
does not apply the SU(3) basis, thus a detailed compar- 
ison with the standard shell model is not easy. 



8 



IV. SUMMARY AND PERSPECTIVES 

Nuclear clustering escapes the description in terms of 
the standard shell model which simply fails to predict 
cluster states at observed low excitation energies around 
cluster-decay threshold. Generic explanation of this enig- 
matic phenomenon involves understanding of the role 
played by the coupling of shell model states via the de- 
cay channel(s). In this work, we have proposed a simple 
approach based on the improved AQCM parametrization 
to describe a transition from the a-cluster wave function 
to the jj-coupling shell model wave function. As com- 
pared to the previous AQCM fl9l - |22j which allows the 
transformation of a single a cluster into a quasi-a clus- 
ter, in the improved AQCM we expand the model space 
and spin directions can be oriented in any direction. In 
the new formulation, all a clusters of a nucleus can be 
changed into quasi-a clusters. Hence, a destructive in- 
terference of the spin-orbit coupling and clustering can 
be investigated in a single variational wave function. 

We applied this improved AQCM wave function to 
12 C. Based on analytical and numerical studies, we have 
shown that the AQCM wave function in the limit R — > 
and A = 1 corresponds to the (si/2) 4 (??3/2) 8 closed shell 
configuration of 12 C. The proposed parametrization al- 
lows to determine an optimal wave function of 12 C in a 
variational procedure for each chosen effective nucleon- 
nucleon interaction. The optimal AQCM ground state 
of 12 C is an intermediate state between the three-a clus- 
ter state and the shell model state with the p^/ 2 subshell 
closure configuration. From a comparison with the AMD 
model, where all nucleons are treated independently, we 
found that the AQCM result is consistent with the AMD 
result (overlap is about 99%) even though the number of 
degrees of freedom in the AQCM trial wave function is 
significantly fewer than in the AMD. 

The AQCM can be applied to heavier nuclei. In sd- 
shell nuclei, a variational wave function has to include 
both clusters and quasi-clusters to describe pure shell 
model configurations. In some cases, the AQCM wave 
function contains even more than one configuration, e.g. 
in 28 Si; one configuration is a pentagon of quasi-a clus- 
ters on the xy-plane and two a clusters along the z-axis, 
and another configuration is a tetrahedron of a clusters 
whose center of gravity is at the origin of the coordi- 
nate system and a triangle of three quasi-a clusters on 
the xy-plane surrounding it. If we take the zero limit 
for the relative distances among a clusters and quasi 
clusters, these two representations become identical and 
both configurations give the lowest shell model configu- 
ration of (s 1/2 ) 4 (p 3/2 ) 8 ( Pl/2 ) 4 (d 5/2 ) 12 at A = 1. On the 
other hand, with finite inter-cluster distances, there is 
difference between these two representations, and these 
two configurations become independent basis states. In 
the low-lying energy region of 28 Si, another configura- 
tion, a 16 0+ 12 C cluster configuration, which has a pro- 



late shape, is also considered to play an important role. 
For such configuration, we can easily prepare the 12 C 
cluster, which is not three a's but the jj-coupling shell 
model one (si/ 2 ) 4 (p3/ 2 ) 8 , using the improved AQCM as 
discussed in the present paper. Utilizing the common 
nucleon-nuclcon interaction, diagonalizing the Hamilto- 
nian matrix comprised of these difference configurations 
is straightforwardly performed using AQCM. 

The AQCM studies open new perspectives in studies 
of cluster-shell competition in various particle-stable sys- 
tems, including neutron-rich carbon nuclei, heavy Na 
nuclei, or hypcrnuclci. The intricate role of dynamical 
symmetries in this phenomenon remains an exciting open 
question for future studies. 

In terms of the shell model approach, the generic expla- 
nation of nuclear clustering, its universal occurrence and 
properties, involves understanding of the role played by 
the coupling of shell model states via the decay chan- 
nels). This coupling leads to the formation of the 
collective near-threshold eigenstate of the open quan- 
tum system. The coupling results in the anti-Hermitian 
component, and interplay between Hcrmitian and anti- 
Hermitian terms is a source of collective near-threshold 
phenomena in the ensemble of shell model states d, Q . 
The coupling to decay channels in AQCM, which could 
be included through the complex scaling method, is a 
challenging subject to be addressed in near future. 
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Appendix A: R expansion of the improved AQCM 
wave function 

In this Appendix, we show that the improved AQCM 
wave function (in this case three quasi-a clusters) be- 
comes the jj-coupling wave function (in this case p 3 / 2 
subshell closure configuration) at the limit of i? — > and 
A = 1. We expand the AQCM wave function in R and 
prove that the i? 8 -order term of the wave function is the 
(si/2) 4 (P3/2) 8 configuration. 

Since proton and neutron parts of the wave function 
are identical, we consider here only the proton part. The 
proton part of the total wave function of the improved 
AQCM with A = 1 is the antisymmctrized product of six 
single particle wave functions: 
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exp 







-v{ ri -R{l,i,0)Y 



,cxp 



cxp 



cxp 



1 J 1 ^ 

- v Y*- R [-2^-— 



(^/ 2 2 )' exp 



where 

(x, y, z) = xe x + ye y + ze z 
stands for the Gaussian center parameters, and 

a A = |t> + 6U) 
is the spin part of the wave function. Here, 



exp 



-v{r x - R(l,i,0)Y 



(A2) 



(A3) 



(A4) 



-v{r 2 -R{l,-i,Q)Y 



-vU-rI- 1 -, 



-V3/2 \ 
1/2 )' 

-1/2 J 



is the spin-up proton, and 



exp 



-v{r 2 -R(l,-i,0)¥ 



, (Al) 



(A5) 



is the spin-down proton of the first quasi-a cluster. Single 
particle wave functions of other four protons are gener- 
ated by rotating both the spatial and spin parts of these 
two orbits about the y-axis by 27r/3 and Att/3. 

Using Eqs. (fT5)l. and ([2"2]l. we expand the proton 
part of the total wave function, Eq. (|A1[) , in R as 



( 2 i ] 



01/2I5, 5 } + a3/2jR| l l ] + 0(jR2) ' a - 1/2 ^' ~5 } + a - 3 / 2Rl l'-l ) + ° (i?2) ' 



ai/2 



III I\ + ^3|I 

2 1 2 ' 2 7 2 '2' 2' 



/ V3.1 1. 1.1 1. 
a -i/2 - — I 2' 2 } + 2 [ 2~2 I 



ai/2 



-1/2 



-i 1 ->+^i---> 

2 1 2 ' 2' 2 '2' 2 7 y 

_V3,1 1, _ 1.1 
2 ' 2' 2 2 2' 2' 



( 1 3 3 3 3 1 3^ 3 1 3V3 3 3.' 

\ / 3^,3 3. 3^3,3 1 3.3 1 1.3 3 

) + a ^ 2R {- — l 2' 2 ] + — l 2> 2 ] ~ 8 1 2' ~2 ] + 8 l r ~2 

( 1.3 3 3 3 1 3^3,3 1 3^3.3 3 \ 

+« 3/2 i? 2> + gig' 2 } - -rt' - ^ + J 

3^,3 3, 3V^,3 1, 3,3 1. 1,3 3 



a-3/2-R 



'2' 2 1 



8 l 2'2 / 



2' 2' 8 1 2 ' 2 



0(i? 2 ), 

+ 0(i? 2 ), 
0(i? 2 ), 
+ 0(R 2 ) 



2v\ 2 999 



7T / 512 



ai/2a-i/2a 3 / 2 a _3/2 



I I\ ,1 _In ,3 1) ,3 I) ,3 _In ,3 _3s 
2 ' 2 /' l 2' 2 7 ' 1 2' 2^' 1 2' 2 ; ' 1 2' 2 ;,l 2' 2' 



0(R 5 ). 



(A6) 
(A7) 



r 



The (si/2) 2 (P3/2) 4 configuration appears as the lowest neutron and proton parts of the wave function have the 
order term in i?, the i? 4 -order term. Because of the anti- same form, the total wave function has the (si/2) 4 (P3/2) 8 
symmctrization, terms up to the i? 3 -order vanish. Since configuration at the i? 8 -order term. 
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